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Abstract 

Using the complex stereographic variable representation for the macrospin, from 
a study of the nonlinear dynamics underlying the generalized Landau-Lifshitz(LL) 
equation with Gilbert damping, we show that the spin-transfer torque is effectively 
equivalent to an applied magnetic field. We study the macrospin switching on a 
Stoner particle due to spin-transfer torque on application of a spin polarized cur- 
rent. We find that the switching due to spin-transfer torque is a more effective 
alternative to switching by an applied external field in the presence of damping. We 
demonstrate numerically that a spin-polarized current in the form of a short pulse 
can be effectively employed to achieve the desired macro-spin switching. 
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1 Introduction 



In recent times the phenomenon of spin-transfer torque has gained much at- 
tention in nanoscale ferromagnets[l,2,3]. Electromigration refers to the recoil 
linear momentum imparted on the atoms of a metal or semiconductor as a large 
current is conducted across. Analogously, if the current is spin-polarized, the 
transfer of a strong current across results in a transfer of spin angular momen- 
tum to the atoms. This has lead to the possibility of current induced switch- 
ing of magnetization in nanoscale ferromagnets. With the success of GMR, 
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this has immense application potential in magnetic recording devices such as 
MRAMs[3,4,5,6]. The phenomenon has been studied in several nanomagnetic 
pile geometries. The typical set up consists of a nanowire[3,7,8,9,10,ll], or a 
spin-valve pillar, consisting of two ferromagnetic layers, one a long ferromag- 
netic pinned layer, and another small ferromagnetic layer or film, separated 
by a spacer conductor layer (see Figure 1). The pinned layer acts as a reser- 
voir for spin polarized current which on passing through the conductor and 
on to the thin ferromagnetic layer induces an effective torque on the spin 
magnetization in the thin film ferro magnet. A number of experiments have 
been conducted on this geometry and the phenomenon has been convincingly 
confirmed [12,13,14,15]. Although the microscopic quantum theory of the phe- 
nomenon is fairly well understood, interestingly the behavior of the average 
spin magnetization vector can be described at the semi-classical level by the 
LL equation with an additional term[16]. 

i z 



film Conductor 

Fig. 1. A schematic diagram of the spin- valve pillar. A thin film ferromagnetic layer 
with magnetization S is separated from long ferromagnetic layer by a conductor. S p 
is the direction of magnetization in the pinned region, which also acts as a reservoir 
for spin polarized current. 

From a different point of view, several studies have focused on magnetic pulse 
induced switching of the macro-magnetization vector in a thin nanodot un- 
der different circumstances [17,18,19,20]. Several experimental studies have 
also focussed on spin-current induced switching in the presence of a magnetic 
field, switching behavior for different choices of the angle of the applied field, 
variation in the switching time, etc., [12,21,22,23,24,25]. A numerical study 
on the switching phenomenon induced by a spin current in the presence of 
a magnetic field pulse has also been investigated very recently in [26]. As an 
extension to two dimensional spin configurations, the switching behavior on a 
vortex has been studied in [27]. 

In this article, by investigating the nonlinear dynamics underlying the gener- 
alized Landau-Lifshitz equation with Gilbert damping, we look at the exciting 
possibility of designing solid state memory devices at the nanoscale, wherein 
memory switching is induced using a spin polarized current alone, without 
the reliance on an external magnetic field. We compare earlier studied switch- 
ing behavior for the macro- magnetization vector in a Stoner particle [17] in 
the presence of an external magnetic field, and the analogous case wherein 
the applied field is now replaced by a spin polarized current induced spin- 
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transfer torque, i.e., with the thin film in the first case replaced by a spin 
valve pillar. It will be shown that a pulse of spin polarized current is more 
effective in producing a switching compared to an applied field. In doing so we 
rewrite the system in terms of a complex stereographic variable instead of the 
macro-magnetization vector. This brings a significant clarity in understanding 
the nonlinear dynamics underlying the macrospin system. Namely, it will be 
shown that, in the complex system, the spin-transfer torque is effectively an 
imaginary applied magnetic field. Thus the spin-transfer term can accomplish 
the dual task of precession of the magnetization vector and dissipation. 

The paper is organized as follows: In Section 2 we discuss briefly the model 
system and the associated extended LL equation. In Section 3 we introduce 
the stereographic mapping of the constant spin magnetization vector to a 
complex variable, and show that the spin-transfer torque is effectively an 
imaginary applied magnetic field. In Section 4 we present results from our 
numerical study on spin-transfer torque induced switching phenomenon of the 
macro-magnetization vector, for a Stoner particle. In particular, we study two 
different geometries for the free layer, namely, (a) an isotropic sphere and (b) 
an infinite thin film. In applications to magnetic recording devices, the typ- 
ical read/write time period is of the order of a few nano seconds. We show 
that, in order to achieve complete switching in these scales, the spin-transfer 
torque induced by a short pulse of sufficient magnitude can be affirmatively 
employed. We conclude in Section 5 with a discussion of the results and their 
practical importance. 



2 The extended LL equation 



The typical set up of the spin-valve pillar consists of a long ferromagnetic 
element, or wire, with magnetization vector pinned in a direction indicated by 
S p , as shown in Figure 1. It also refers to the direction of spin polarization of 
the spin current. A free conduction layer separates the pinned element from 
the thin ferromagnetic film, or nanodot, whose average spin magnetization 
vector S(t) (of constant magnitude So) is the dynamical quantity of interest. 
The cross sectional dimension of the layers range around 70 — lOOnm, while 
the thickness of the conduction layer is roughly 2 — 7nra[3,20]. The free layer 
thus acts as the memory unit, separated from the pinned layer cum reservoir 
by the thin conduction layer. It is well established that the dynamics of the 
magnetization vector S in the film in the semiclassical limit is efficiently de- 
scribed by an extended LL equation[16]. If m(= {mi,m 2 ,m 3 } = S/Sq) is the 
unit vector in the direction of S, then 

dm -± ^ dm . ^ ^ - , ^ . ^ A . . . 

— = -7m x H eff + Am x — - -fag(P, m ■ S p )m x (m x S p ), (1) 
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a = 2SjTe (2) 

Here, 7 is the gyromagnetic ratio (= 0.0176 Oe^ns^ 1 ) and So is the satura- 
tion magnetization (Henceforth we shall assume 47iSo = 8400, the saturation 
magnetization value for permalloy). The second term in (1) is the phenomeno- 
logical dissipation term due to Gilbert [28] with damping coefficient A. The last 
term is the extension to the LL equation effecting the spin-transfer torque, 
where A is the area of cross section, j is the current density, and V is the 
volume of the pinned layer, 'a', as defined in (2), has the dimension of Oe, and 
is proportional to the current density j. g(P,m- Sp) is given by 

^ * ' ^ = /(P)(3 + l.S p )-4 ; fiP) = (3> 

where f(P) is the polarization factor introduced by Slonczewski [1], and P(0 < 
P < 1) is the degree of polarization of the pinned ferromagnetic layer. For 
simplicity, we take this factor g to be a constant throughout, and equal to 1. 
H e ff is the effective field acting on the spin vector due to exchange interaction, 
anisotropy, demagnetization and applied fields: 

Heff ^exchange ~\~ Hanisotropy ~\~ H demagnetization ~\~ H applied-, (4) 



where 

^■exchange = -DV 2 lh, (5) 

H anisotropy = k{u! • e||)e||, (6) 

V • H demagnetization = ~ 47rS'oV • ih. (7) 

Here, k is the strength of the anisotropy field, ey refers to the direction of 
(uniaxial) anisotropy, In what follows we shall only consider homogeneous 
spin states on the ferromagnetic film. This leaves the exchange interaction 
term in (4) redundant, or D = 0, while (7) for H demagnetization is readily solved 
to give 

— * 

H demagnetization = -4irS (A^iTOiX + N 2 m 2 y + N 3 m 3 Z,), (8) 

where JVj, i = 1,2,3 are constants with Ni + N 2 + N 3 = 1, and {x, y, z} are the 
ortho normal unit vectors. Equation (1) now reduces to a dynamical equation 
for a representative macro-magnetization vector m. 



In this article we shall be concerned with switching behavior in the film purely 
induced by the spin-transfer torque term, and compare the results with earlier 
studies on switching due to an applied field [17] in the presence of dissipation. 
Consequently, it will be assumed that H ap piied = in our analysis. 
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3 Complex representation using stereographic variable 



It proves illuminating to rewrite (1) using the complex stereographic variable 
tt defined as[29,30] 

n 

so that 

n + n 
mi = TTW ; m2 

For the spin valve system, the direction of polarization of the spin-polarized 
current S p remains a constant. Without loss of generality, we chose this to be 
the direction z in the internal spin space, i.e., S p = z. As mentioned in Sec. 
2, we disregard the exchange term. However, for the purpose of illustration, 
we choose H app i ie( i = {0, 0, h a3 } for the moment but take h a3 = in the later 
sections. Defining 

e|| = {sin#j| cos0||, sin 6\\ sin^y, cos fly} (11) 
and upon using (9) in (1), we get 



uii + im 2 
1 + m 3 



—i 



. (n - n) 
'i + |fi| 2 



m 3 



i + \n\ 2 ' 



(9) 
(10) 



(1 -i\)Cl= -7(0 - ih a3 )Q + im\\K-f\ cos 0\\Q - ]- sin 0y (e^n 



fi 2 e-^i) 



(i+W) 



N 3 {1 - \Q\ 2 )Q - — (1-Q 2 
2 



_^(1 + Q2_| Q |2 )Q _ 



(iVi 



- I^l 2 )fi 



(12) 



where my = m • ey. Using (10) and (11), my, and thus (12), can be written 
entirely in terms of Q. 

It is interesting to note that in this representation the spin-transfer torque 
(proportional to the parameter a) appears only in the first term in the right 
hand side of (12) as an addition to the applied magnetic field h a3 but with a 
prefactor —i. Thus the spin polarization term can be considered as an effective 
applied magnetic field. Letting k — 0, and iVi = iV 2 = N 3 in (12), we have 



(1 — i\)Cl = — 7(0 — ih a3 )fl, 
which on integration leads to the solution 



(13) 



ft(t) = fi(0) exp (-(a - ih a3 )~ft/(l - 

= Q(0) exp(- 1 + , A | 2 70 exp (-2 - - 7 t). (14) 



5 



The first exponent in (14) describes relaxation, or switching, while the second 
term describes precession. From the first exponent in (14), we note that the 
time scale of switching is given by l/(a+Xh a3 ). A being small, this implies that 
the spin-torque term is more effective in switching the magnetization vector. 
Further, letting h a3 = 0, we note that in the presence of the damping term 
the spin transfer produces the dual effect of precession and dissipation. 

To start with we shall analyze the fixed points of the system for the two cases 
which we shall be concerned with in this article: (i) the isotropic sphere char- 
acterized by Ni = N 2 = N 3 = 1/3, and (ii) an infinite thin film characterized 
by Nx = = N 3 , N 2 = 1. 

(i) First we consider the case when the anisotropy field is absent, or k = 0. 
From (12) we have 



(1 - = - ai Q - n4 ffi 2 [N 3 (l - \Vl\ 2 )Q -^(1-Q 2 - \Q\ 2 )Q 

W-"»>fl].(i5) 

In the absence of anisotropy (k = 0), we see from (15) that the only fixed 
point is Qq — 0. To investigate the stability of this fixed point we expand (15) 
up to a linear order in perturbation 5Q around f2 . This gives 

(l-i\)6(l = -a75fi-i747r5o[iV3-i(iV 1 +iV 2 )]5fi+i727r 1 So(iV 1 -A^ 2 )(5n. (16) 

For the isotropic sphere, Ni — N 2 — N 3 — 1/3, (16) reduces to 

(1 - iX)5tl = -a^Sn. (17) 

We find the fixed point is stable since a > 0. For the thin film, JVi = = 
N 3 ,N 2 = I. (16) reduces to 

(1 - iX)5Q = -a-fSn + i~f2irS 5n - i-firS 5tt. (18) 

This may be written as a matrix equation for \l/ = (Sfl , 5Q) T , 

^ = M^, (19) 

where M is a matrix obtained from (18) and its complex conjugate, whose 
determinant and trace are 

II 1 + A 2 ' V ) 1 + A 2 V ) 

Since |M| is positive, the fixed point Q = is stable if TV|M| < 0, or, 
(a + 27r,SoA) > 0. 
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The equilibrium point (a), Qo = 0, corresponds to rh = z. Indeed this holds 
true even in the presence of an applied field, though we have little to discuss 
on that scenario here. 

(ii) Next we consider the system with a nonzero anisotropy field in the z 
direction. (12) reduces to 



[I - = -a-fQ + iK'y- ~ 



2^ ... 

nm - \n\ 2 )n 



(i + |fi| 2 ) (i + |fi| 2 ) 

_ _ j n |2 )n _ ^ (1 + n 2 _ | n) 2 )n _ (^1-^ 0] . (21) 

Here again the only fixed point is f2 = 0. As in (i), the stability of the 
fixed point is studied by expanding (21) about fl to linear order. Following 
the same methodology in (i) we find the criteria for stability of the fixed 
point for the isotropic sphere is (a + Xk) > 0, while for the thin film it is 
(a + \(K + 2nS )) > 0. 

(iii) With nonzero k in an arbitrary direction the fixed point in general moves 
away from z. 

Finally, it is also of interest to note that a sufficiently large current leads to 
spin wave instabilities induced through spin-transfer torque [31,32]. In the 
present investigation, however, we have assumed homogeneous magnetization 
over the free layer, thus ruling out such spin wave instabilities. Recently we 
have investigated spin wave instabilities of the Suhl type induced by an applied 
alternating field in thin film geometries using stereographic representation [30]. 
It will be interesting to investigate the role of a spin-torque on such instabil- 
ities in the spin valve geometry using this formulation. This will be pursued 
separately. 



4 Spin-transfer torque induced switching 



We now look at the interesting possibility of effecting complete switching of the 
magnetization using spin-transfer torque induced by a spin current. Numerical 
studies on switching effected on a Stoner particle by an applied magnetic field, 
or in the presence of both a spin-current and applied field, in the presence of 
dissipation and axial anisotropy have been carried out recently and switching 
has been demonstrated [17,26]. However, the intention here is to induce the 
same using currents rather than the applied external fields. Also, achieving 
such localized magnetic fields has its technological challenges. Spin-transfer 
torque proves to be an ideal alternative to accomplish this task since, as we 
have pointed out above, it can be considered as an effective (albeit complex) 
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magnetic field. In analogy with ref. [17], where switching behavior due to an 
applied magnetic field has been studied we investigate here switching behav- 
ior purely due to spin-transfer torque, on a Stoner particle. Numerical results 
in what follows have been obtained by directly simulating (12) and making 
use of the relations in (10), for appropriate choice of parameters. It should 
be remembered that (12) is equivalent to (1), and so the numerical results 
have been further confirmed by directly numerically integrating (1) also for 
the corresponding parameter values. We consider below two samples differing 
in their shape anisotropies, reflected in the values of (Ni, N 2 , N 3 ) in the de- 
magnetization field: a) isotropic sphere, Ni = N 2 = N 3 = 1/3 and b) a thin 
film Ni = = N$, N 2 = 1. The spin polarization S p of the current is taken 
to be in the z direction. The initial orientation of m is taken to be close to 
— z. In what follows this is taken as 170° from z in the (z — x) plane. The 
orientation of uniaxial anisotropy ey is also taken to be the initial direction 
of m. With these specified directions for S p and ey the stable fixed point is 
slightly away from z, the direction where the magnetization m is expected to 
switch in time. A small damping is assumed, with A = 0.008. The magnitude 
of anisotropy k is taken to be 45 Oe. As stated earlier, for simplicity we have 
considered the magnetization to be homogeneous. 



4-1 Isotropic sphere 



It is instructive to start by investigating the isotropic sphere, which is char- 
acterized by the demagnetization field with Ni = N 2 = N 3 = 1/3. With these 
values for (N 1: N 2 , jV 3 ), (12) reduces to 



(1 -i\)tl = -a7Q + im||re7[cos0||ft- ^ sin 6»|| (e^H -Q 2 e^)] (22) 

A constant current of a = 10 Oe is assumed. Using (2), for typical dimen- 
sions, this equals a current density of the order 10 8 A/cm 2 . We notice that 
for the isotropic sample the demagnetization field does not play any role in 
the dynamics of the magnetization vector. In the absence of anisotropy and 
damping the spin-transfer torque term leads to a rapid switching of S to the 
z direction. This is evident from (22), which becomes 

SI = -a 7 fi, (23) 

with the solution = f2 e"" 7 *, and the time scale for switching is given by 
1/07. Figure 2. a shows the trajectory traced out by the magnetization vector 
S, for 5 ns, initially close to the — z direction, switching to the z direction. 
Figure 2.b depicts the dynamics with anisotropy but no damping, all other 
parameters remaining same. While the same switching is achieved, this is 
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more smoother due to the accompanying precessional motion. Note that with 
nonzero anisotropy, z is not the fixed point any more. The dynamics with 
damping but no anisotropy (Figure 2.c) resembles Figure 2. a, while Figure 
2.d shows the dynamics with both anisotropy and damping. 





Fig. 2. Trajectory of the magnetization vector m, obtained by simulating (12) for 
the isotropic sphere (N\ = N2 = N3 = 1/3), and using the relations in (10), for 
a = 10 Oe (a) without anisotropy and damping, (b) with anisotropy but no damping, 
(c) without anisotropy but nonzero damping and (d) with both anisotropy and 
damping nonzero. The results have also been confirmed by numerically integrating 
(1). The arrows point in the initial orientation (close to — z) and the direction of the 
spin current z. Evolution shown is for a period of 5ns. Note that the final orientation 
is not exactly z in the case of nonzero anisotropy ((b) and (d)). 



It may be noticed that Figures 2.c and 2.d resemble qualitatively Figures 
2. a and 2.b, respectively, while differing mainly in the time taken for the 
switching. It is also noticed that switching in the absence of anisotropy is 
faster. Precession assisted switching has been the favored recording process in 
magnetic memory devices, as it helps in keeping the exchange interaction at 
a minimum [18, 19]. The sudden switching noticed in the absence of anisotropy 
essentially refers to a momentary collapse of order in the magnetic media. 
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This can possibly lead to strong exchange energy and a breakdown of our 
assumption regarding homogeneity of the magnetization field. However, such 
rapid quenching assisted by short high intensity magnetic pulses has in fact 
been achieved experimentally [33]. 

A comparison with reference [17] is in order. There it was noted that with an 
applied magnetic field, instead of a spin torque, a precession assisted switch- 
ing was possible only in the presence of a damping term. In Section 3 we 
pointed out how the spin transfer torque achieves both precession and damp- 
ing. Consequently, all four scenarios depicted in Figure 2 show switching of 
the magnetization vector without any applied magnetic field. 

4-2 Infinite thin film 

Next we consider an infinite thin film, whose demagnetization field is given by 
Ni = = N 3 and N 2 = 1. With these values (12) becomes 



Here again in the absence of anisotropy Q = is the only fixed point. Thus the 
spin vector switches to z in the absence of damping and anisotropy (Figure 
3a). In order to achieve this in a time scale of 5ns, we find that the value of a 
has to be of order 50 Oe. Again the behavior is in stark contrast to the case 
induced purely by an applied field [17], wherein the spin vector traces out a 
distorted precessional trajectory. As in Sec. 4.1, the trajectory traced out in 
the presence of damping is similar to that without damping (Figure 3c). The 
corresponding trajectories traced out in the presence of anisotropy are shown 
in Figures 3b and 3d. 

4-3 Switching of magnetization under a pulsed spin-polarized current 

We noticed that in the absence of uniaxial anisotropy, the constant spin polar- 
ized current can effect the desired switching to the orientation of S p (Figure 
2). This is indeed the fixed point for the system (with no anisotropy). Fig- 
ure 2 traces the dynamics of the magnetization vector in a period of 5 ns, in 
the presence of a constant spin-polarized current. However, for applications in 
magnetic media we choose a spin-polarized current pulse of the form shown 
in Figure 4. It may be recalled here that, as was observed in 4.2, with a spin 



(1 



i\)Ct = — cryfi + im\\Kj cos9\\Q 



2 



sin0||(e^ii -fi 2 e-^ii) 




(24) 
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Fig. 3. Trajectory of the magnetization vector S in a period of 5ns, obtained as 
earlier by numerically simulating (12), and also confirming with (1), with demagne- 
tization field with by N\ = = N^, N2 = 1 and a = 50 Oe, and all the parameter 
values are as earlier. As in Figure 1, the S and initial orientation are indicated 
by arrows, (a) Without anisotropy or damping, (b) nonzero anisotropy but zero 
damping, (c) without anisotropy but nonzero damping and (d) both anisotropy and 
damping nonzero. As earlier, in the presence of nonvanishing anisotropy, the fixed 
point is not the z axis. 

polarized current of sufficient magnitude, the switching time can indeed be 
reduced. We choose a pulse, polarized as earlier along the z direction, with 
rise time and fall time of 1.5ns, and a pulse width, denned as the time interval 
between half maximum, of 4ns. We assume the rise and fall phase of the pulse 
to be of a sinusoidal form, though, except for the smoothness, the switching 
phenomenon is independent of the exact form of the rise or fall phase. 



In Figures 5 and 6, we show trajectories of the spin vector for a period of 
25 ns, for the two different geometries, the isotropic sphere and a thin film. 
The action of the spin torque pulse, as in Figure 4, is confined to the first 
5 ns. We notice that, with the chosen value of a, this time period is enough 
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1 2 3 4 5 

time (ns) 

Fig. 4. Pulse form showing the magnitude of a, or effectively the spin-polarized 
current. The rise and fall phase are assumed to be of a sinusoidal form. The rise 
and fall time are taken as 1.5ns, and pulse width 4 ns. The maximum magnitude 
of a is 150 Oe. 



to effect the switching. In the absence of anisotropy, the direction of S p is 
the fixed point. Thus a pulse of sufficient magnitude can effect a switching in 
the desired time scale of 5ns. From our numerical study we find that in order 
for this to happen, the value of a has to be of order 150 Oe, or, from (2), 
a current density of order 10 9 Aj cm 2 , a magnitude achievable experimentally 
(see for example [34]). Comparing with sections 4.1 and 4.2, we note that the 
extra one order of magnitude in current density is required due to the duration 
of the rise and fall phases of the pulse in Figure(4). Here again we contrast 
the trajectories with those induced by an applied magnetic field [17], where 
the switching could be achieved only in the presence of a uniaxial anisotropy. 



In Figure 5b for the isotropic sphere with nonzero crystal field anisotropy, 
we notice that the spin vector switches to the fixed point near z axis in the 
first 5 ns. However the magnetization vector precesses around z after the 
pulse has been turned off. This is because in the absence of the spin-torque 
term, the fixed point is along en, the direction of uniaxial anisotropy. Due to 
the nonzero damping term, the spin vector relaxes to the direction of ey as 
time progresses. The same behavior is noticed in Figure 6b for the thin film, 
although the precessional trajectory is a highly distorted one due to the shape 
anisotropy. 
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Fig. 5. Evolution of the magnetization vector S in a period of 25ns induced by the 
spin-polarized current pulse in Figure 4, (a) with and (b) without anisotropy for 
the isotropic sample, with N\ = N2 = A3 = 1/3 all other parameters remaining 
same. A nonzero damping is assumed in both cases. The current pulse acts on the 
magnetization vector for the first 5 ns. In both cases switching happens in the first 
5 ns. In the presence of nonzero anisotropy field, (b), the magnetization vector 
precesses to the fixed point near z. 




Fig. 6. Evolution of the magnetization vector S in a period of 25ns induced by 
the spin-polarized current pulse in Figure 4, for a infinite thin film sample, with 
Ni = = A3, and N2 = 1, all other parameters remaining same, along with a 
nonzero damping, (a) Without anisotropy and (b) with anisotropy. As in Figure 5, 
switching happens in the first 5 ns. 

5 Discussion and conclusion 



We have shown using analytical study and numerical analysis of the nonlinear 
dynamics underlying the magnetization behavior in spin-valve pillars that a 
very effective switching of macro-magnetization vector can be achieved by a 
spin transfer-torque, modeled using an extended LL equation. Rewriting the 
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extended LL equation using the complex stereographic variable, we find the 
spin-transfer torque term indeed acts as an imaginary applied field term, and 
can lead to both precession and dissipation. It has also been pointed out why 
the spin-torque term is more effective in switching the magnetization vector 
compared to the applied field. On application of a spin-polarized current the 
average magnetization vector in the free layer was shown to switch to the 
direction of polarization of the spin polarized current. For a constant current, 
the required current density was found to be of the order of 10 8 Aj cm 2 . For 
recording in magnetic media, switching is achieved using a stronger polarized 
current pulse of order 10 9 A/ cm 2 . Currents of these magnitudes have been 
achieved experimentally. 
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